We identify two distinct low-energy sectors in the classi- 
cal isotropic antiferromagnetic Heisenberg spin-S chain. In 
the continuum limit, we show that two types of rotation gen- 
erators arise for the field in each sector. Using these, the 
Lagrangian for sector I is shown to be that of the nonlinear 
sigma model. Sector II has a null Lagrangian; Its Hamilto- 
nian density is just the Pontryagin term. Exact solutions are 
found in the form of magnons and precessing pulses in I and 
moving kinks in II. The kink has 'spin' S. Sector I has a higher 
minimum energy than II. 
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There are certain interacting many-body systems, which 
display a close correspondence between their extreme quan- 
tum behavior and strictly classical behavior, in the following 
sense: If the quantum energy spectrum of the one-dimensional 
lattice Hamiltonian is exactly solvable, then the classical con- 
tinuum limit of its dynamical equation, which is usually a non- 
linear partial differential equation, is completely integrable, 
supporting soliton solutions. As examples M we note the 
correspondence between : the quantum spin-^ isotropic fer- 
romagnetic chain and the classical Landau-Lifshitz equation; 
bosons with delta-function interactions and the nonlinear 
Schrodinger equation. 

Like its ferromagnetic (FM) counterpart in the example 
cited above, the isotropic antiferromagnetic (AFM) Heisen- 
berg chain described by the Hamiltonian 

H = —J ) ] S B ' Sn+l] (1) 
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where S„ denotes the spin at the n th site, (S n ) 2 = S 2 and 
J < 0, is also exactly solvable M in the extreme quantum case 
S = |. Faddeev and Takhtajan || used Bethe ansatz tech- 
niques to show that the low-lying excitation here is a quan- 
tum kink with spin 1. Subsequently, using a formulation 
proposed by Mikeska Haldane |J showed that the clas- 
sical continuum dynamics of this system can be mapped to 
that of the 0(3) nonlinear sigma model (NLSM) Lagrangian, 
whose equations of motion are known to be integrable jf|. 
However, as he rightly pointed out, a finite anisotropy term 
must be added to the AFM Hamiltonian in (1), to obtain a 
classical kink solution (i.e., a topological soliton) for the spin 
configuration. (As is obvious from the functional form of this 
solution, it ceases to be a kink in the isotropic limit.) Soon af- 
terwards, the quantum, large spin S, AFM chain was mapped 
[0 , H in the continuum to the Lagrangian of the NLSM with 
a topological term. Since this term is a 'total divergence', 
the equations of motion are unaffected by it; and therefore, 
here too, the presence of an anisotropy appears to be crucial 
to support a kink. Given the quantum<->classical correspon- 
dence mentioned in the beginning, and the rigorously-proven 
existence of the quantum kink in the isotropic case, this ap- 
parent absence of the classical kink in the isotropic limit is 
strange indeed This motivates us to take a deeper look 
into the classical equations of motion generated by the basic 
AFM Hamiltonian (1). 

In this Letter, in addition to showing how a kink arises 
in the isotropic AFM as an exact dynamical solution in the 
continuum, we present several other new results that provide 
important insights into the low-energy dynamical structure 
underlying this system: On the discrete lattice, we identify 
two distinct low-energy sectors I and II that one must ana- 



lyze for a full understanding of this system. We derive the 
continuum dynamical equations of (1), taking care to ensure 
that the Poisson bracket structure is preserved exactly. These 
equations help us to directly identify the corresponding total 
angular momentum in each sector. A careful scrutiny of the 
discrete system shows that in the continuum limit, there are 
in fact two types of rotation generators that get associated 
with the vector field (in each sector). When appropriately 
taken into account, these yield the above total angular mo- 
mentum. Using this, its Lagrangian is shown to be that of 
the NLSM. In this sector, we obtain exact nonlinear magnons 
and precessing pulse solutions for the unit vector field. Sector 
II is described by a single equation for the vector field. The 
associated Hamiltonian turns out to be a homogeneous func- 
tion of the first degree in the canonical momenta. Using a 
Lagrangian analog of Dirac's [jioj formulation of constrained 
dynamics, we show that the Lagrangian of this sector vanishes 
identically. We find exact moving kinks for the field equation. 
This kink does not require an anisotropy to support it, in con- 
trast to the NLSM sector I. We show that the total angular 
momentum of this kink is simply S. It carries a 'vorticity' 
with it, reminiscent of the superfluid vortex. Typically, the 
minimum energy of sector I is greater than that of sector II. 
We begin our analysis by defining 

n(2i + 1/2) = [S(2i + 1) - S(2i)]/2S; (2) 

l(2i + 1/2) = [S(2i + 1) + S(2z)]/2a. (3) 

Here a is the lattice separation. Eqs. (2) and (3) are just the 
classical counterparts of the 'staggered magnetization' n, and 
the 'spin density' 1, customarily used in the quantum formu- 
lation [^|. The dynamical equations for S(2i) and S(2i + 1) 
are found from dS n /dt = {S n ,H}, where {, } denotes the 
Poisson bracket (PB). 5*" satisfy the usual angular momen- 
tum algebra {S%, St} = S nm e a/3j SZ. Using this, Eq.(l) 
yields dS n /dt — JS n x (S n+ i + S n _i). In an AFM chain, 
the low-energy dynamics corresponds to [S(2i + 1) + S(2i)] 
being very small. In the continuum limit a — * 0, this leads 
to two possibilities : (I) [S(2i + 1) + S(2i)] = o(a m ),m > 1, 
and (II) [S(2i + 1) + S(2i)] = 0, identically. Noting that the 
usual AFM spin wave velocity c — 2| J\Sa must remain finite 
as a — > 0, our detailed analysis |nj shows that in (I), only 
m — 1 leads to nontrivial dynamics. All m > 2 lead to trivial 
spin configurations, constant in space and time. (II) corre- 
sponds to a 'dimer-like' configuration in which every second 
pair of spins on the chain gets 'locked' in exactly antiparal- 
lel directions, and has interesting dynamics. We proceed to 
analyze these two cases. 

Sector (I): Here, l(2i + i) is nonzero. The two constraints 
n 2 = 1 + (a/S) 2 l 2 and n ■ 1 = ensure that there are four 
independent variables per unit cell (2a), describing this sector 
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as required (l2|. In the continuum limit a — > 0, we have, 
n 2 — > 1, (2i + 5)0 — > x, 8 m n/2a — > 5(a; — y), etc. Using the 
angular momentum algebra of the spins S„, one obtains for 
the continuum, 

{l a (x),f(y)} = e a ^F(x)S(x-y) ; {n a {x), f{y)} = 

= e a *n'(x)8(x-v) ; {n a (x), n (y)} -> 0. (4) 

Thus the spin density 1 behaves like an angular momentum 
for the unit vector field n in the continuum limit. From the 
discrete equations of motion for S(2i + 1) and S(2i) arising 
from (1), we write down the equations for l(2i + i) and n(2i + 
|), and find their continuum limit by using appropriate Taylor 
expansions for l(a;±2a) and n(x±2a). We obtain the following 
coupled equations : 

n t = (2 C /S)[(l-| n!C )xn] (5) 

l t = c[(l - Sn*) x n]„. (6) 

The subscripts x and t stand for partial derivatives. Eq. (5) 
immediately yields 

l=(S/2)[i(nxn,)+n.] (7) 

c 

On substituting Eq. (7) in Eq. (6), we get 

n x [(nti/c 2 ) - n xx ] = 0. (8) 

Eq. (8) is just the Euler-Lagrange equation for the NLSM, 
and is integrable ||. Using Eqs. (2) and (3) in (1), and taking 
the continuum limit, 

H= (c/S) J ' dx[{\~^n x ) 2 + (S 2 /4)n 2 )] = J h dx. (9) 

To find the Lagrangian density C for (9), we need to identify 
the total canonical momentum 7T. Now, an inspection of Eq. 
(5), which essentially is an equation of motion for a 'nonlinear 
rigid rotator' n, hints strongly that its total angular momen- 
tum is in fact Ai = (1 — |ni) rather than just 1, which satisfies 
(4). To understand the origin of this fact, we first define 

L(2i+i) = [S(2i-l) + S(2z) + S(2i + l) + S(2z + 2)]/4a. (10) 

This is an obvious extension of 1 defined in (3) and repre- 
sents the spin density at (2i + |) due to a cluster of 4 spins 
around this point, with the appropriate denominator 4a. A 
short calculation using the usual PB algebra obeyed by spins 
S n shows that in the continuum limit, L also behaves like 
an angular momentum of n, i.e., its components satisfy (4). 
Next, writing the expressions for the various spins appearing 
in (10) in terms of n and 1 by using (2) and (3), and taking 
the continuum limit, we get 

L = l-Sn x (11) 

Further, let us denote 1 in (3)) by L2, and L in (10) by L4, with 
the subscripts denoting the number of spins in the cluster, and 
similarly define spin densities L6,Ls,Lio etc., corresponding 
to clusters of 6, 8, 10. .spins. Then, computing the PBs of these 
L n with n(2i + §) shows that in the continuum limit, each 
of them satisfies Eq. (4) , and so can qualify as the angular 



momentum of n. In addition, a tedious calculation shows that 
in the continuum, the spin densities corresponding to Lg, L10, 
L14 etc, reduce simply to 1, while L4, Ls, L12 etc, reduce to L 
given in (11). Thus the contributions from the successive L n 's 
alternate between 1 and L, leading to an average spin density 
(1 + L)/2. On using Eq. (11), this indeed gets identified with 
the the total angular momentum Ax: 

Ai = (1 + L)/2 = (1 - Jn,) = (S/2c)(n x n t ), (12) 

where we have used Eq. (7) to write the last equality. 

To find the Lagrangian, n is written in spherical polar co- 
ordinates as 

n = (sin 6 cos <j>, sin 9 sin <j>, cos 9) (13) 

Let 7Tfl and ir^, denote the canonical momenta satisfying 
{9(x),TT 9 (y)} = -0(x),7i>(y)} = S(x - y). Noting that 
(n, lie, n^fsin 9)~ 1 )) form a unit orthogonal triad, we have 
as usual M, Ai = (n x [none + 7r^n^(sin6) -2 )]. Here, the 
subscripts 9 and (j> of n denote partial derivatives. Com- 
paring this with the last equation in (12), we get -Kg = 
{S/2c)9 t and ti> = (S/2c) sin 2 6<j> t . Further, from Eq. (9), 
h = [ng 2 + n (t> 2 (sm9)~ 2 + sin 2 9(f> 2 x + 9% Using these in £ 
=TVg9t + it<p4>t — h, we obtain £=[(^-)n 2 — n 2 ]. This is just the 
NLSM Lagrangian. Note that both the rotation generators 1 
and L must be used for the correct identification of the to- 
tal angular momentum of n in the continuum, which in turn 
is needed to find the associated Lagrangian. This is a sub- 
tle point, indeed. Our detailed analysis shows that ignoring 
either one of these and treating the other alone as the total 
angular momentum would lead to a spurious topological term 
in the classical NLSM Lagrangian. 

Analyzing Eq. (8) in terms of 9 and <f> variables using 
Eq. (13), we find exact traveling wave solutions of the form 
9 — 0(ujit — kix) and (f> = {tut — kx) . Details jll| will be 
published elsewhere. For these solutions, the product of the 
velocities of the 9-wave and the </>-wave turns out to be vgv,/, 
=(u>/fc)(wi/fci)=c 2 , implying that these waves travel in the 
same direction. As an example, we have the following 'pre- 
cessing pulse' that travels with velocity c: 

n = (—tanh^ cos£,— tanh£ sin£, sech£), (14) 

where £ = (cut — kx) and ui — ck. As |x| — > 00, although 
— > 0, indicating a 'pulse', n' 1 ' and rS 2 ^ do not vanish, 
but precess on a plane. We also find exact nonlinear analogs 
of the usual small amplitude spin wave solutions of Eq. (8). 
Typically, the traveling waves in this sector have a constant 
energy density. The consequences of this will be discussed at 
the end. 

Sector (II): As already mentioned, here we have S(2i) = 
-S(2i + 1). Thus Eq.(2) yields: 

n(2i + 1/2) -> S(2i + 1)/S= -S(2i)/S = N(2i + 1/2) (15) 

and Eq.(3) gives 1(2? + |) = 0. From Eq.(15), it is easily ver- 
ified that in the continuum limit a — > 0, {N a (x), N°(x)} — > 
2§-5{x ~ y)e al3 ~< N"< (x) -> 0. Thus N(x) behaves like a unit 
vector field. Using the dynamical equations for S(2i) and 
S(2i+1) leads in the continuum limit to the following single 
equation for N: 

N( = c (N x N x ) (16) 
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Since N 2 = 1, writing it in terms of polars (6,(j>) as we did 
for n in Eq. (13), leads to: 

t = —csw.0<p x ■ <j> t = c0 x /sm.9 (17) 

Using Eq.(15) in Eq.(l) yields the following continuum Hamil- 
tonian: 

= J h dx (18) 

Proceeding as in sector I leads to the total angular momentum 
A 2 = -(S/2)N X for this sector. Further, A 2 = (N x [pgNg + 
p</,N^(sin6>)~ 2 )]. Writing N, = Ng9 x + N<j,<fi x , we compare 
the above two expressions of A2 to obtain the momenta 

pe = -(5/2) sin 9 4> x ■ P4 , = (3/2) sin 9 6 X . (19) 

Note the dependence of the momenta on the 'generalized coor- 
dinates', 9 and <j> and their spatial derivatives. It is clear that 
to generate the equations of motion (17), from 9 t = SHo/Spg 
and cj>t = SHo/Sptf,, the corresponding Hamiltonian density ho 
must have the form 

ho = -c [p e sin 9 <j> x - (p^/ sin 9) 9 X ] (20) 

Thus ho is homogeneous of the first degree in the p's. Note 
that there are no independent equations for (pe)t and (p</,)t in 
this sector. A short calculation shows that the corresponding 
Lagrangian density Co=Pe9t +P(p4>t — ho vanishes identically. 
This may be surprising at first sight, but not if it is realized 
that the dynamics in this sector is just an example of the La- 
grangian analog of Dirac's ^) constrained dynamics, encoun- 
tered in particle mechanics. Treating the dynamical equations 
(16) as constraints, which leads to N 2 /c 2 = N 2 , our null 
Lagrangian may be written in the form £<p = (N 2 /c 2 — N 2 ), 
which yields pe — (S/2c)0t and p$ = (5/2c) sin 2 0(f>t. These 
yield Eq.(19), on using Eq. (17). Further, using these in (20) 
gives ho = (S/2)sm9[9 x </> t - 0t<t> x \ = (5/2)N • (N, x N t ). 
Thus the energy density in this sector is effectively the Pon- 
tryagin (topological) term. To our knowledge, this is the first 
example in condensed matter physics where Dirac's method 
for constrained dynamics finds an explicit application. 

We have found jllj exact traveling wave solutions of the 
form 9 — 9(ui\t — kix) and (f> = (uit — kx) for the dynamical 
equations (17) in this sector. In general these are traveling 
kinks, for which the product vgv^ — (ui/k)(uii/ki)——c 2 , im- 
plying that these two waves travel in opposite directions to 
each other, in contrast to sector I. As an example, we present 
the following solution: 

N= (sech^i cos£,sec/i£i sin£, — tanh^i). (21) 

Here £1 = (cot + kx), f = (uit — kx) and ui = ck. Since 
N = S(2i + l)/5 = -S(2i)/5, Eq. (21) represents a kink 
profile in each sublattice, moving with velocity c, with the 
spins precessing around the z-axis with velocity — c. The odd 
site spins rotate from 9 — at x — > — 00 to 9 = n at x — > +00, 
and vice versa for the even site spins. Unlike in sector I, the 
kink is a topological soliton that interpolates between two dis- 
tinct Neel states. Typically, the traveling wave solutions in 



this sector have an energy density which is a 'lump', vanishing 
as I a; I — * 00. The total energy of the kink is found from Eq. 
(18) to be E K = 2Sui = 2Sck = 2Sc/F, where V = (1/k) is 
the width of the kink. Thus the energy vanishes as k — > 0. 
The z-component of the integrated total angular momentum 
of the kink is M {3) = / A { 2 3) dx = (-5/2) / N x 3) dx = 5, while 
Mjp , i — 1,2 vanish. Thus the total angular momentum of 
the kink is 'self-generated' and can hence be interpreted as 
its 'spin'. Its magnitude is simply the spin S of the chain (1). 
Further, there is a clear spatial separation between the 'vortic- 
ity' field A2 ( - i \x,t) = 4£ sech 2 (kx + uit) which is a pulse, and 
the kink field N^(x,t) = —tanh(kx + uit) , as in a superfluid. 

The above results correspond to [S(2i) + S(2i + 1)] = 0, 
implying l(2i + i) = (see Eq.(3)). If we consider the other 
possibility [S(2i) + S(2i - 1)] = 0, then clearly L(2i + ±) = 
(see Eq.(10)). Repeating the same steps as before, we obtain 
= — c(N x Nz), This implies waves traveling with velocity 
c opposite in direction to that in Eq.(16 ). Further, 1 = SN X 
here, yielding A 2 = (5/2)Na,. Thus Sector II supports kink 
solutions with velocity ± c depending on which of the 'dimer' 
bonds is in the antiparallel-locked configuration. These kinks 
have identical energies, total angular momenta etc. 

Returning to Sector I, we have shown that the isotropic 
AFM chain maps to the NLSM model without a topolog- 
ical term. This appears to be in agreement with the re- 
sult that Haldane jB| had obtained in the isotropic limit, 
in his seminal paper using Mikeska's ansatz M. It is 
therefore instructive to check if this ansatz, namely, S„ = 
(— l)"5'(sin& Tl cos(^„,sin6 n sin(^ n ,cos6 n ), with 9„ — 9(x) + 
a(—l) n a(x), 4> n = 4>(x) + a(—l) n f5(x), x — na satisfies the 
basic relationship 1 = | [in x nt + n x ] (see Eq.(7)) derived 
from the exact continuum equations of (1)). A long calcula- 
tion shows that a and /3 must equal (9 x /2) and (<j> x /2) re- 
spectively, to satisfy the above relationship. This in turn can 
be shown to imply that the ansatz is indeed consistent for 
1=0 (Sector II) and therefore leads to our Eq. (17), which 
supports kink solitons. Using this equation, the Lagrangian 
(Eq.(10) in [B]) is easily seen to vanish in the isotropic limit, 
as expected. 

Eq. (17) also supports instanton solutions, since it is just 
the (1+1)D analog of the 2D Belavin-Polyakov equation [|13j. 
This equation has been shown to be integrable by mapping it 
to the elliptic Liouville equation |l4|. Thus the NLSM equa- 
tion (Eq.(8)) of sector I and the Belavin-Polyakov equation 
(Eq.(16)) of sector II are both integrable. 

Finally, we discuss the energy associated with the two sec- 
tors: For sector I, a short calculation using Eq. (9) shows that 
the minimum energy density corresponds to l m in = (5/2)^ 
yielding h m in = (c/S)!^^. Recalling the restriction that 
both generators of rotation, 1 and L have to be finite in this 
sector, leads to the result that h m in is finite in this sector. In 
contrast, for sector II, from Eq. (18), the minimum energy 
density is directly seen to be ho,min = 0, corresponding to 
= 0, implying that it is continuously connected to the 
ground state energy. Thus sector I has a higher minimum 
energy than sector II. 

All our results so far have been fully classical. Although 
a rigorous quantum mechanical analysis of the two sectors 
is needed for the full understanding of the energy spectrum, 
semiclassically, we see that for sector I, our classical result 
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implies E m in ~ 1(1 + 1). Since we have the restriction I ^ 0, 
only the next higher integer value of I is allowed, leading to a 
gap. On the other hand, there is no such restriction on either 
of the angular momenta in sector II, and hence E m in ~ 
here. Thus at any finite temperature, the ground state of the 
isotropic AFM chain gets disordered due to the excitation of 
kinks rather than magnons ||. Further, as we showed earlier, 
the spin of the kink is S. Specializing to the S = | case, our 
analysis suggests that a kink in the isotropic AFM chain is 
a spin-| entity which is gapless, in agreement with Faddeev 
and Takhtajan's result || as well. 
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